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Abstract 

A  theoretical,  two-dimensional,  along-the-channel  model  has  been  developed  to  design  fuel  channels  for  proton  exchange  membrane 
(PEM)  fuel  cells.  This  has  been  implemented  by  solving  the  resultant  ordinary  differential  equation  with  a  straightforward  shooting 
computational  scheme.  With  such  a  design  tool,  an  analysis  can  be  made  of  the  effects  due  to  some  operation  and  design  parameters, 
such  as  inlet  velocity,  inlet  pressure,  catalyst  activity,  height  of  channel,  and  porosity  of  gas-diffusion  layer  to  obtain  a  fuel  cell  with  high 
performance.  Present  results  indicate  that  there  is  always  a  trade-off  between  higher  power  density  and  higher  efficiency  of  the  fuel  cell. 
Namely,  a  design  for  higher  power  density  (a  better  performance)  is  always  accompanied  with  a  higher  fuel  efficiency  (or  a  larger  fuel 
consumption  rate  and  a  higher  fuel  cost),  and  vice  versa.  When  some  relevant  physical  parameters  are  determined  experimentally  and 
applied  in  the  present  model,  a  quantitative  design  for  a  fuel  cell  of  either  high  efficiency  or  high  performance  is  feasible. 

©  2003  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  proton  exchange  membrane  (PEM)  fuel  cell  is  an 
electrochemical  device  which  combines  fuel  (hydrogen)  and 
oxidant  (oxygen)  to  produce  electricity,  and  water  and  heat 
are  the  major  by-products.  In  past  decades,  substantial  ef¬ 
forts  [1-4]  have  been  devoted  to  reducing  the  cost  as  well  as 
promoting  the  efficiency  of  the  fuel  cell.  In  this  respect,  the 
design  of  high-efficiency  fuel  channels  is  one  of  the  impor¬ 
tant  issues  [5-19].  Several  different  types  of  fuel  channels 
have  been  used  in  practical  designs,  such  as  straight  chan¬ 
nels,  surpentine  channels,  and  interdigitate  channels.  The 
morphology  of  the  channel  also  varies,  namely,  meander, 
spiral  or  straight  types.  For  a  meander  or  a  spiral  channel, 
the  length  may  be  several  meters  long.  In  such  a  long  chan¬ 
nel,  the  fuel  may  be  entirely  consumed  before  exit,  which 
implies  that  a  certain  portion  of  the  cell  may  not  have  fuel 
for  chemical  reactions.  That  is,  a  part  of  the  cell  may  not 
produce  any  electrons  during  operation  and  this  reduces  the 
efficiency.  Therefore,  analysis  of  the  channel  flow  becomes 
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a  necessity  in  fuel  cell  design  and  an  efficient  and  conve¬ 
nient  theoretical  model  for  channel  analysis  is  essential. 

To  date,  there  have  been  two  major  approaches  for 
the  analysis  of  channel  flow.  One  uses  computational 
fluid  dynamic  (CFD)  techniques  to  examine  the  two-  or 
three-dimensional  flow  in  fuel  channels  [5-10],  the  other 
uses  a  one-dimensional  approximation  approach  to  inves¬ 
tigate  the  variation  of  flow  structure  and  includes  the  fuel 
concentration  as  well  as  the  current  generated  along  the 
channel  [11-19].  In  the  CFD  approaches,  mathematical 
models  are  usually  developed  for  the  whole  PEM  fuel  cell, 
which  may  consider  the  conservation  of  mass,  momen¬ 
tum  and  energy,  equations  governing  the  electrochemical 
reaction,  and  various  kinds  of  physical  properties  of  the 
components  such  as  diffusivity  of  the  gas-diffusion  layer, 
electro-osmosis  in  membrane,  fuel  convection  across  the 
membrane,  chemical  reaction  and  activity  in  catalyst  lay¬ 
ers,  and  membrane  hydration  [20-23].  With  this  approach, 
however,  the  computation  is  rather  time-consuming  and  the 
analysis  procedure  is  so  tedious  that  the  computation  of 
the  whole  flow  field  in  channels  becomes  inefficient  and  is 
sometimes  even  an  unnecessary  step  to  obtain  the  fuel  cell 
design. 
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Nomenclature 

c  dimensionless  concentration 

C  concentration  (mol  m-3) 

Co  inlet  concentration  (molm-3) 

F  Faraday  constant  (96,500  C  mol- 1 ) 
h  channel  height  (m) 

i  local  current  density  (Am-2) 

10  inlet  current  density  (A  m-2) 

Ki  i  =  1-5,  constants  in  Eq.  (15) 

L  channel  length  (m) 

M  molar  weight  (kg mol-3) 

n  transferred  electron  number 

P  pressure  (Pa) 

P0  inlet  pressure  (atm) 

P  pressure  gradient  (Pam-1) 

R  universal  gas  constant  (8.314  Jmol-1  K-1) 

T  gas  temperature  (K) 

u  velocity  in  v-direction  (ms-1) 

uo  inlet  velocity  in  v-direction  (ms-1) 
u  cross-sectional  averaged  velocity  in  v-direction 

u  dimensionless  velocity  in  v-direction 

v  velocity  in  y-direction  (ms-1) 

Vo  suction  velocity  (ms-1) 

Greek  letters 

a  charge  transfer  coefficient 

y  reaction  order 

0  over  potential  (V) 

11  viscosity  (kgm-1  s-1) 

§  empirical  constant  regarding  the  slip  condition 
at  the  porous  boundary 
p  density  (kgm-3) 


A  more  efficient  and  convenient  scheme,  which  can  to 
some  extent  reach  the  goal  of  fuel  cell  design  both  quantita¬ 
tively  and  qualitatively,  is  the  channel  model  [11-19].  Since 
the  flow  passes  the  fuel  channel  very  rapidly,  cross-sectional 
variations  of  flow  structure  and  other  physical  parameters 
such  as  fuel  concentration,  fuel  density  and  fuel  tempera¬ 
ture  can  be  ignored.  A  major  concern  is  the  variation  of 
the  relevant  physical  properties  along  the  channel.  It  was 
found  that  the  along-the-channel  model  is  much  simplified 
and  that  the  equations  can  be  solved  more  conveniently  and 
efficiently.  Recently,  several  along-the-channel  models  have 
been  developed  for  the  above-mentioned  purpose.  Nguyen 
and  co-workers  [10,11]  have  proposed  a  set  of  governing 
equations  which  include  water  and  energy  transport  across 
the  membrane  and  have  considered  heat  removal  along  the 
channel.  Other  workers  [12-14]  have  examined  the  gas  dy¬ 
namics,  concentration  decay  and  current  drop  along  the 
channel.  Argyropoulos  et  al.  [15-19]  have  investigated  the 
pressure  drop  and  the  temperature  variation  along  the  chan¬ 
nel  by  considering  the  mass  and  energy  conservation  of 


two-phase  flows.  In  these  studies,  however,  the  channel  flow 
was  assumed  to  have  a  constant  velocity  along  the  channel. 
Under  most  circumstances,  however,  this  assumption  may 
result  in  significant  discrepancies  with  reality. 

In  the  present  work,  a  two-dimensional  theoretical  model 
is  developed,  which  includes  a  continuity  equation,  momen¬ 
tum  equations  and  the  Tafel  equation.  At  the  bottom  of  the 
channel,  the  consumption  of  fuel  due  to  chemical  reactions 
is  modeled  by  way  of  sucking  the  fuel  through  the  porous 
boundary.  This  model  allows  investigation  of  the  variations 
of  the  flow  structure,  the  fuel  concentration,  and  the  current 
density  along  the  channel.  A  systematic  parametric  study 
is  implemented  to  examine  the  influence  of  relevant  design 
and  operation  parameters  on  the  effective  length  of  the  fuel 
channel.  Specifically,  the  effects  due  to  the  inlet  velocity 
and  pressure  and  the  porosity  of  the  gas-diffusion  layer  at 
the  bottom  can  be  investigated.  The  activity  of  the  catalyst 
attached  to  the  gas  diffusion  layer  can  also  be  analyzed. 


2.  The  theoretical  model 


A  two-dimensional  horizontal  channel  flow  is  shown 
schematically  in  Fig.  1.  The  channel  has  a  constant  height 
and  is  sufficiently  long  that  the  conditions  at  both  the  en¬ 
trance  and  the  exit  do  not  affect  the  flow  of  the  domain  under 
consideration.  The  fluid,  either  hydrogen  in  the  anode  or 
oxygen  in  the  cathode,  is  assumed  to  be  an  ideal  gas.  Since 
the  flow  velocity  is  so  high,  the  variation  of  temperature 
along  the  channel  is  assumed  to  be  negligible.  It  is  also  as¬ 
sumed  that  the  bipolar  plate  in  which  the  fuel  channel  is  built 
is  an  ideal  collecting  electrode  with  no  ohmic  loss.  The  over¬ 
potential  can  then  be  maintained  at  a  constant  value  along  the 
channel  [12].  As  a  result,  the  local  current  density  is  a  func¬ 
tion  of  the  fuel  concentration  at  each  specific  position  [7]. 
For  such  a  chemically  active  flow,  the  continuity  equation  is: 


dp  d(pu)  d(pv) 

d  t  dx  +  dy 


in  which  u  =  u(x,  y)  and  v  =  v(x,  y)  are  the  velocity  in 
the  v-and  y-directions,  respectively.  Due  to  the  fact  that 
the  channel  height  is  very  small  compared  with  the  chan¬ 
nel  length,  it  may  be  assumed  that  the  density  distribution 


Fig.  1.  Schematic  description  of  flow  considered.  Two-dimensional  viscous 
flow  is  bounded  between  bipolar  plate  (solid  wall)  at  top  and  gas-diffusion 
layer  (porous  wall)  at  bottom.  Velocity  profile  is  essentially  a  parabolic 
curve.  Because  of  the  porous  wall  there  is  a  velocity  slip  at  bottom. 
Because  of  chemical  reaction  in  the  gas-diffusion  layer,  there  is  a  suction 
flow  across  the  porous  wall,  accounted  for  by  suction  velocity  Vq(x). 
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across  the  channel  is  uniform,  i.e.,  p  =  p(x ).  Because 
p(x)  =  MC(x ),  Eq.  (1)  can  be  converted  to: 


3(0/)  d(Cv) 
dx  +  3  y 


Integrating  Eq.  (2)  along  y  and  applying  Leibnitz’s  rule 
gives: 


3  rh  dc  rh  rhd v 

C—  udy+  —  l  udy  -\-  C  I  —  dy  =  0. 
dxj o  dx  Jo  Jo  9y 


By  letting  u(x)  =  (l//z) JqU(x,  y)dy  be  the  averaged  ve¬ 
locity  over  the  cross-section  of  the  channel  and  by  applying 
the  non-slip  condition  at  the  top  u(x ,  h)  =  0  and  the  slip 
condition  at  the  bottom  u(x ,  0)  =  § u(x ),  we  obtain: 


h 


d[C(x)  u(x)] 
dx 


—  C(x)  v(x,  0)  =  —C(x)  Vo(x), 


(4) 


where  Vo  is  the  suction  velocity  along  the  bottom  of  chan¬ 
nel,  which  is  defined  as  Vo  =  (1  /nF)(i(x)/C(x)).  Thus, 
Eq.  (4)  becomes: 

d \C(x)u(x)]  1  i(x) 


dx 


h  nF 


(5) 


The  relation  between  the  overpotential  rj  and  the  electric 
current  density  i  is  governed  by  the  Tafel  equation: 


i(x)  =  io 


C(x) 

Co 


y 


exp 


aF 

RT 


T] 


(6) 


so  that  the  suction  velocity  Vo  can  be  expressed  by: 


Vo(x)  = 


Note  that,  the  slip  boundary  condition  at  the  bot¬ 
tom  u(x,  0)  =  %u(x)  is  derived  from  the  so-called 

Beavers-Joseph  boundary  condition  at  the  interface  be¬ 
tween  a  fluid  and  a  porous  layer  [24],  and  the  parameter 
§  essentially  accounts  for  the  permeability  of  the  porous 
media  since  the  factor  due  to  the  velocity  gradient  at  the 
bottom  is  absorbed  into  the  averaged  velocity  u(x).  Given 
that  the  velocity  gradient  does  not  change  significantly 
along  the  channel,  a  change  in  §  may  be  seen  as  a  change  in 
the  permeability  or,  equivalently,  a  change  in  the  porosity 
of  the  porous  medium  below.  In  the  present  study,  it  con¬ 
sidered  that  §  varies  from  0.1  to  1  [24]  due  to  the  fact  that 
the  porosity  of  the  gas-diffusion  layer  at  the  bottom  is  high. 

For  the  momentum  equation,  it  is  assumed  that,  for  the 
present  two-dimensional  channel  flow,  the  velocity  in  the 
v-direction  is  much  larger  than  that  in  the  y-direction.  After 
applying  order  analysis  on  the  momentum  equations  (or  the 
Navier-Stokes  equations)  in  both  v-and  y-directions,  a  single 
momentum  equation  results  as  follows: 


p(x) 


3  u  3  u 

u - 1-  v  — 

dx  dy 


d  P 
dx 


The  pressure  gradient  along  the  channel  is  assumed  to  be 
constant,  i.e.,  (dP/dx)  =  constant  =  P.  Under  the  normal 


operating  conditions  of  a  1  kW  fuel  cell,  P  =  12 Pam-1 
[2,12-14],  which  is  to  be  used  in  the  present  analysis.  Inte¬ 
grating  Eq.  (8)  along  the  height  yields: 


ph 

3  u  3  u 

du 

MCI 

u - 1-  v  — 

dy  =  —hP  +  ll  — 

Jo 

dx  dy 

dy 

To  simplify  Eq.  (9)  further,  it  is  assumed  that  u(x ,  y)  is  a 
quasi-parabolic  velocity  profile  [25]  defined  as  u(x,  y)  = 
A(x)y2  +  B(x)y  +  G(x).  At  y  =  0,  u(x,  0)  =  ^u(x)  and 
this  leads  to  G(x)  =  at  y  =  h,  u(x,  y)  =  0  and 

B(x)  =  — ha(x )  —  (%/  h)u(x).  To  obtain  A(v),  the  non- slip 
boundary  condition  is  applied  at  the  top  u(x,  h)  =0  and  this 
results  in  A(x)  =  —  (6/  /z2)  ( 1  —  (1/2 )%)u(x).  Consequently 
the  approximated  velocity  function  in  v-direction  is  obtained 
as: 


u(x,  y)  =  u(x) 


2 ■ 
h 

(10) 


This  equation  implies  that  the  horizontal  velocity  is  a 
second-order  parabolic  function  of  y,  with  a  small  slip  at 
the  bottom  of  channel,  as  shown  schematically  in  Fig.  1. 
To  obtain  v(x,  h),  Eq.  (10)  is  substituted  in  Eq.  (8)  and  the 
resultant  equation  is  integrated  along  y  and  the  boundary 
conditions  of  v  are  applied  at  the  top  and  the  bottom  of  the 
channel,  i.e.,  v(x,  h)  =  0  and  v(x,  0)  =  Vq(x),  yielding: 


v(x,y)  = 


[C(x)u(x)  f 
C(x) 
-V0(x). 


r 


h 2 


(ii) 


Eqs.  (10)  and  (11)  are  substituted  in  Eq.  (9)  to  give  the 
following  equation  for  the  channel  flow: 


2  ?  1  6 
—f  -  -H  +  - 

15s  5s  5 


+  —  +  12  1 
M 


d[C(x)u2(x)]  % 


dx 


+  -C(x)  u(x)  Vo (x) 
h 


2fy  Mh 2 


/x 

u(x)  ~  0. 


(12) 


Eqs.  (5),  (6)  and  (12)  are  the  governing  equations  for  the 
two-dimensional  flow  along  the  channel,  in  which  the  vari¬ 
ations  of  the  fuel  concentration  C(x),  the  velocity  u(x)  and 
the  current  density  i(x)  are  to  be  solved.  Since  these  equa¬ 
tions  have  an  initial  value  problem,  the  initial  velocity  at  the 
entrance,  uq,  and  the  initial  concentration  fed  to  the  chan¬ 
nel,  Co,  are  required.  It  is  therefore  assumed  that  u(x)  = 
u(x)/uo  and  c(x)  =  C(x)/Co  and  these  two  relationships 
are  substituted  into  the  equations  to  give: 


d(cu) 

dx 


+  K\cy  =  0, 


(13) 


dx 


T  H-  K/[Cyu  T  K$u  —  0, 


(14) 
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where 

_  ip  (aF  \ 

1  nFhC0u0  CXP  yHT11)' 

K2  =  fe2-&  +  i 

_  p 

—  MCqUq  ’  (15) 

z'o  / otF  \ 

K4  ~  ^ nFhC0u0  CXP  yRT11)  ’ 

K5  =  n(\-  -A  — Jf — . 

V  2  7  Mh2C0UQ 

The  initial  conditions  become  u  =  1  and  c  =  1 .  These  two 
equations  are  solved  by  a  fourth  order  Runge-Kutta  scheme. 
Eqs.  (13)  and  (14)  are  the  two  simplified  equations  for  the 
averaged  horizontal  velocity  u  and  averaged  fuel  concentra¬ 
tion  c.  After  obtaining  these  two  values,  the  Tafel  equation 
(Eq.  (6))  can  be  applied  to  obtain  the  local  current  den¬ 
sity  along  the  channel.  The  combination  of  these  equations 
and  the  shooting  scheme  becomes  a  convenient  tool  to  de¬ 
termine  various  physical  parameters  relevant  to  the  design 
of  fuel  channels  of  PEM  fuel  cells.  This  scheme  is  differ¬ 
ent  from  those  developed  in  other  studies  [11-19]  in  which 
the  velocity  variations  in  both  x-  and  y-directions  were  ig¬ 
nored  and  only  the  effective  length  was  considered.  Using 
this  simplified  model  allows  not  only  the  examination  of  the 
effective  length  of  the  fuel  channel  under  various  operation 
conditions,  but  also  the  effects  due  to  relevant  design  and 
operation  parameters  on  the  fuel  cell  performance. 

3.  Typical  flow  field:  an  example 

Analyses  in  the  present  and  the  following  sections  are 
made  on  the  basis  of  the  base  case;  the  values  of  its  relevant 
physical  parameters  are  shown  in  Table  1.  This  case  essen¬ 
tially  corresponds  to  a  1  kW  PEM  fuel  cell  under  normal 
operation  conditions  [2,12-14].  The  results  in  terms  of  the 
variations  of  velocity,  fuel  concentration  and  current  den¬ 
sity  along  the  channel  are  presented  in  Fig.  2.  It  is  found 
that  the  fuel  concentration  decays  monotonically  along  the 
channel  due  to  the  chemical  reaction  occurring  at  the  bot¬ 
tom  of  channel,  while  the  flow  velocity  increases  along  the 
channel  because  of  the  depletion  of  fuel  downstream.  Both 
effects  are  the  direct  consequence  of  fuel  consumption  along 
the  channel,  which  in  turn  leads  to  a  decrease  in  generated 
current  (Fig.  2(b))  because  of,  again,  the  depletion  of  fuel 
downstream.  These  findings  suggest  that  variations  of  these 
physical  properties  occur  simultaneously,  and  therefore  will 
be  considered  as  a  whole  instead  of  separately  as  in  previ¬ 
ous  studies  [11-19]  Also  the  velocity  is  taken  as  constant 
while  the  fuel  concentration  changes  along  the  channel. 

The  generation  of  current  density  decreases  along  the 
channel  because  of  the  decay  in  fuel  concentration  (see 
Fig.  2(b)).  Note  that,  for  the  present  base  case,  the  fuel 


Table  1 

Values  of  physical  parameters  of  base  case  corresponding  to  1  kW  PEM 
fuel  cell  [2,12-14]. 


Channel  side 

Cathode 

Gas  flow 

Oxygen 

Half-reaction 

02(g)  +  4H+  +  4e“  - 

2H20 

Channel  temperature,  T  (K) 

353.15 

Inlet  gas  velocity,  uq  (ms-1) 

0.1 

Inlet  gas  pressure,  Po  (atm) 

2 

Inlet  gas  concentration,  Co  (molm-3) 

69.00 

Exchange  current  density,  4)  (Am-2) 

10“5 

Activation  overpotential,  ij  (V) 

0.3 

Reaction  order,  y 

0.5 

Electrons  transferred  in  reaction,  n 

4 

Charge-transfer  coefficient,  a 

2.0 

Molar  weight,  M  (kg  mol-1) 

32  x  10“3 

Viscosity,  /x  (kg  ms-1) 

2  x  10“5 

Channel  height,  h  (m) 

10“3 

Channel  length,  L  (m) 

0.1 

Pressure  gradient,  P  (Pam-1) 

10 

Universal  gas  constant,  R  (Jmol-1  K-1) 

8.314 

Faraday  constant,  F  (Cmol-1) 

96500 

Slip  velocity  fraction,  £ 

0.1 

K\  (m-1) 

13.77 

k2 

1.181 

K3  (m-1) 

452.9 

K4  (m-1) 

1.377 

K5  (m-1) 

1032.6 

concentration  decrease  by  about  50%  and  the  local  current 
density  by  about  25%  within  the  10  cm  long  channel.  The 
variations  in  these  two  parameters  are  very  significant  for 
such  a  short  distance,  implying  that  an  accurate  calculation 
of  the  effective  channel  length  cannot  be  ignored  in  the 
design  of  a  high-efficiency  fuel  cell.  In  a  practical  sense, 
too  long  a  channel  will  result  in  a  large  dead  region  in  the 


Fig.  2.  Along-the-channel  variations  of  major  quantities  of  present  prob¬ 
lem:  (a)  dimensionless  velocity  u  and  concentration  ~c\  (b)  generated  cur¬ 
rent  density.  Results  are  calculated  on  basis  of  base  case  shown  in  Table  1 . 
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bipolar  plate,  while  too  short  a  channel  will  cause  a  great 
part  of  the  fuel  to  leave  the  cell  without  reaction.  A  long 
channel  gives  a  low  local  current  density  downstream  (or  a 
low  performance  of  the  fuel  cell)  but,  with  respect  of  fuel 
consumption,  a  long  channel  length  may  ensure  that  the 
fuel  is  consumed  before  leaving  the  exit.  The  latter  will 
result  in  a  better  efficiency  of  fuel  usage. 

4.  Physical  parameter  effects 

Several  physical  parameters  have  significant  effects  on  the 
channel  flow,  which  may  in  turn  affect  both  the  design  and 
operation  conditions  of  the  fuel  cell.  These  physical  param¬ 


eters  include  the  flow  velocity  at  the  inlet,  the  fuel  concen¬ 
tration,  the  activity  of  the  catalyst,  the  channel  height,  and 
the  porosity  of  the  gas-diffusion  layer.  The  following  exam¬ 
ines  the  physical  effects  due  to  these  five  parameters  with 
special  attention  to  their  influence  on  the  effective  length  of 
the  channel.  When  the  effect  of  a  particular  physical  param¬ 
eter  is  considered,  this  parameter  will  be  changed  systemat¬ 
ically  while  the  other  four  physical  parameters  are  fixed  at 
the  values  shown  in  Table  1. 

4.1.  Effect  of  inlet  velocity 

The  variation  of  flow  velocity,  fuel  concentration  and  local 
current  density  along  the  channel  are  shown  in  Fig.  3(a)-(c), 


u 


(a) 


(b)  0-50 


(A  cm'2) 


(c) 


Fig.  3.  Effects  of  inlet  velocity  on  along-the-channel  variations  of:  (a)  dimensionless  velocity  u;  (b)  dimensionless  concentration  c;  (c)  local  current  density. 
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respectively,  for  variation  in  inlet  velocity  from  0.07  to 
lms-1.  For  a  lower  inlet  velocity,  because  the  fuel  has  more 
time  to  diffuse  into  the  active  layer  at  the  bottom  to  react, 
the  fuel  concentration  decays  more  rapidly  (Fig.  3(b)).  De¬ 
crease  in  fuel  concentration  leads  to  an  increase  in  flow  ve¬ 
locity  in  accordance  with  the  conservation  of  mass  (Fig.  3(a)) 
and  results  in  a  larger  decrease  in  the  local  current  density 
due  to  the  higher  depletion  of  fuel  or  a  lower  reaction  rate 
(Fig.  3(c)).  It  is  interesting  to  note  that,  for  an  inlet  veloc¬ 
ity  of  1  ms-1,  the  changes  in  flow  velocity,  fuel  concentra¬ 
tion  and  local  current  density  along  the  channel  are  all  very 
small.  As  a  result,  practically,  if  a  uniformly  high  current 


density  along  the  channel  (and  thus  a  higher  power  density) 
is  required,  then  a  larger  inlet  fuel  velocity  can  be  applied 
to  the  cell,  but  at  the  expense  of  a  higher  fuel  consumption 
rate.  On  the  other  hand,  if  the  fuel  efficiency  is  the  major 
concern,  it  is  necessary  to  apply  a  lower  inlet  velocity  so 
that  the  fuel  consumption  (and  thus  the  reaction)  along  the 
channel  can  be  implemented  more  completely. 

4.2.  Effect  of  fuel  concentration 

From  the  equation  of  state  of  an  ideal  gas,  i.e.,  Co  = 
Po/RT ,  the  fuel  concentration  is  a  function  of  the  partial 


Fig.  4.  Effects  of  inlet  fuel  pressure  on  along-the-channel  variations  of:  (a)  dimensionless  velocity  u\  (b)  dimensionless  concentration  c;  (c)  local  current 
density. 
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pressure  of  the  fuel  gas.  Accordingly,  the  change  in  fuel  con¬ 
centration  can  be  attributed  to  the  change  in  partial  pressure 
of  the  fuel.  Thus,  to  examine  the  effects  of  fuel  concentra¬ 
tion,  the  inlet  pressure  Pq  of  the  fuel  was  changed  from  1 .5 
to  5  atm.  The  results  are  presented  in  Fig.  4(a)-(c)  and  show 
that  for  a  higher  inlet  pressure  (or  a  higher  fuel  concentra¬ 
tion),  and  thus  a  larger  pressure  gradient  along  the  channel 
when  the  pressure  at  outlet  is  assumed  to  be  constant,  the 
increase  in  flow  velocity  along  the  channel  is  smaller  (but 
the  overall  velocity  is  higher)  (Fig.  4(a))  because  the  con¬ 
sumption  rate  of  the  fuel  along  the  channel  is  smaller  (and 
the  decrease  of  the  fuel  concentration  along  the  channel  is 
smaller)  (see  Fig.  4(b)).  Due  to  the  smaller  fuel  consumption 


rate,  the  decrease  in  current  density  along  the  channel  is  also 
smaller.  As  a  result,  to  ensure  that  the  current  density  along 
the  channel  can  be  uniformly  high,  or  equivalently  to  have  a 
fuel  cell  of  higher  power  density,  the  cell  must  be  supplied 
with  a  fuel  of  higher  concentration  (a  higher  inlet  pressure) 
but,  again,  at  the  expense  of  a  higher  fuel  consumption  rate. 

4.3.  Effect  of  catalyst  activity 

The  activity  of  the  catalyst  can  be  indicated  by  the  value  of 
the  overpotential  77,  which  can  be  converted  into  the  current 
density  i  by  the  Tafel  equation  (Eq.  (6)).  The  present  analy¬ 
sis  considers  three  different  activities,  corresponding  to  three 
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Fig.  5.  Effects  imposed  by  different  catalysts  on  along-the-channel  variations  of:  (a)  dimensionless  velocity  u;  (b)  dimensionless  concentration  c;  (c) 
local  current  density. 
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overpotentials  r\  =  0.28,  0.30  and  0.32  V.  The  corresponding 
current  densities  are  z’o  =  10-4,  10~5  and  10~6  A  cm-2.  A 
catalyst  of  higher  activity  leads  to  a  faster  reaction  rate  and 
thus  a  faster  speed  of  electron  generation  and,  accordingly 
corresponds  to  a  lower  overpotential  loss  and  a  higher  cur¬ 
rent  density.  Note  that,  a  different  catalyst  activity  implies 
a  different  fuel  consumption  rate,  which  can  be  reflected  by 
the  change  in  suction  velocity  Vo  at  the  bottom  of  the  chan¬ 
nel  (see  Eq.  (7)). 

The  results  shown  in  Fig.  5(a)-(c)  illustrate  that  the  differ¬ 
ence  in  activity  of  the  catalyst  examined  in  the  present  study 
does  not  result  in  an  obvious  difference  in  either  the  flow 
velocity  (Fig.  5(a))  or  the  fuel  concentration  (Fig.  5(b)),  be¬ 
cause  the  curves  of  different  activities  virtually  overlap  each 
other.  On  the  other  hand,  it  has  a  significant  effect  on  the  lo¬ 


cal  current  density,  as  shown  in  Fig.  5(c),  namely:  a  higher 
catalyst  activity  leads  to  a  higher  local  current  density  and 
results  in  a  slightly  more  rapid  decrease  along  the  channel. 
This  indicates  that  an  efficient  scheme  to  raise  the  fuel  cell 
performance  without  consuming  more  fuel  is  to  use  a  cata¬ 
lyst  of  higher  activity.  This  is  a  common  scenario  accepted 
by  fuel  cell  researchers  world- wide.  Accordingly,  the  devel¬ 
opment  of  a  high-activity  catalyst  for  PEM  fuel  cells  has 
been  a  major  research  issue,  or  may  be  the  most  important 
one. 

4.4.  Effect  of  channel  height 

Since  it  is  assumed  that  the  fuel  is  well  mixed  across  the 
channel  height,  a  larger  height  of  channel  is  equivalent  to  a 
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Fig.  6.  Effects  of  channel  height  on  along-the-channel  variations  of:  (a)  dimensionless  velocity  u;  (b)  dimensionless  concentration  c;  (c)  local  current 
density. 
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greater  mass  of  the  fuel  per  unit  length,  or  a  more  sufficient 
supply  of  the  fuel.  Five  different  heights  of  channel  are  con¬ 
sidered,  viz.,  h  =  0.5,  0.75,  1.0,  1.25  and  1.5  mm;  the  results 
are  shown  in  Fig.  6(a)-(c).  It  is  seen  that  a  smaller  height 
of  channel  leads  to  a  larger  increase  in  flow  velocity  and  a 
larger  decrease  in  fuel  concentration  along  the  channel.  Al¬ 
ternatively,  a  smaller  height  of  channel  results  in  a  larger 
decrease  of  local  current  density  along  the  channel  and  thus 
a  smaller  power  density  (or  a  worse  performance)  for  the 
fuel  cell.  Note  that,  since  both  the  flow  velocity  and  the  fuel 
concentration  have  been  normalized  by  the  corresponding 
inlet  conditions,  a  larger  decrease  of  fuel  concentration  in 
fact  accounts  for  a  smaller  amount  of  fuel  consumed  along 
the  channel.  According  to  these  results,  the  application  can 
be  made  as  follows.  To  have  a  fuel  cell  of  higher  fuel  ef¬ 
ficiency,  it  is  necessary  to  have  a  larger  height  of  channel 
because  the  fuel  can  be  consumed  more  efficiently.  Corre¬ 
spondingly,  to  have  a  fuel  cell  of  higher  power  density,  it  is 


necessary  to  have  a  larger  height  of  channel  so  that  a  larger 
amount  of  fuel  can  be  supplied. 

4.5.  Effect  of  porosity  of  gas -diffusion  layer 

The  change  in  porosity  of  the  gas-diffusion  layer  can  be 
reflected  in  the  present  model  by  changing  the  parameter  § 
at  the  bottom  boundary.  Physically,  as  mentioned  above,  a 
larger  §  means  a  larger  porosity  of  the  gas-diffusion  layer. 
Five  cases  are  considered,  namely,  §  =  0.1,  0.2,  0.3,  0.5  and 
1.0.  The  results  show  that  a  larger  §  (or  a  larger  porosity) 
results  in  a  smaller  increase  in  flow  velocity  (Fig.  7(a))  as 
well  as  a  smaller  decrease  in  fuel  concentration  (Fig.  7(b)) 
because  a  larger  portion  of  fuel  is  allowed  to  be  sucked  into 
the  porous  boundary  at  the  bottom  due  to  the  larger  poros¬ 
ity.  Accordingly,  a  greater  chemical  reaction  occurs  in  the 
gas-diffusion  layer  and  leads  to  a  uniformly  high  current  den¬ 
sity  along  the  channel  (Fig.  7(c)).  Accordingly,  as  discussed 


(b) 


Fig.  7.  Effects  of  porosity  of  gas-diffusion  layer  accounted  for  by  parameter  £  on  along-the-channel  variations  of:  (a)  dimensionless  velocity  u\  (b) 
dimensionless  concentration  ~c\  (c)  local  current  density. 
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in  above,  a  fuel  cell  with  higher  power  density  requires  a 
larger  porosity  of  the  gas-diffusion  layer  to  maintain  a  uni¬ 
formly  high  current  density  along  the  channel.  This  conclu¬ 
sion  is  consistent  with  previous  analyses  of  the  optimization 
of  the  porosity  of  the  gas-diffusion  layer  ([26], and  references 
wherein)]  which  show  that  a  higher  porosity  gives  better  per¬ 
formance  of  the  fuel  cell  because  in  the  gas-diffusion  layer 
a  larger  space  (or  pores)  is  given  to  the  water  transportation 
so  that  the  flooding  can  be  prevented  in  high  power  density 
regimes.  Thus,  a  larger  porosity  of  the  gas-diffusion  layer  is 
preferred  for  a  fuel  cell  with  a  high  power  density. 

5.  Concluding  remarks 

A  theoretical  model  for  two-dimensional  flow  in  fuel 
channels  has  been  developed.  This  allows  calculation  of 
variations  in  fuel  velocity,  fuel  concentration,  and  current 
density  along  the  channel.  The  resultant  ordinary  equations 
and  the  initial  conditions  at  the  inlet  of  the  channel  consist 
of  an  initial  value  problem,  which  can  be  easily  solved  by 
a  straightforward  shooting  scheme.  The  combination  of  the 
simplified  equation  model  and  the  popular  shooting  scheme 
becomes  a  convenient,  as  well  as  an  efficient,  scheme  for 
the  design  of  fuel  cell  channels,  such  that  several  design  and 
operation  parameters  can  be  determined  for  a  fuel  cell  with 
either  high  power  density  or  high  fuel  efficiency. 

To  have  a  fuel  cell  of  high  power  density,  results  obtained 
in  the  present  study  suggest  that  it  is  necessary  to:  (1)  in¬ 
crease  the  inlet  velocity  of  fuel;  (ii)  increase  the  inlet  pres¬ 
sure  of  fuel;  (iii)  decrease  the  height  of  fuel  channel;  and/or 
(iv)  increase  the  porosity  of  the  gas-diffusion  layer.  These 
four  schemes  also  apply  to  fuel  cell  with  longer  fuel  chan¬ 
nels  (or  a  larger  size).  Nevertheless,  a  high  power  density  is 
always  followed  with  a  high  fuel  consumption  rate,  which 
results  in  a  low  fuel  efficiency.  In  other  words,  to  enhance 
the  fuel  efficiency  of  a  cell  or  to  apply  a  fuel  cell  with  shorter 
fuel  channels  (or  smaller  size),  it  is  necessary  to  apply  a 
lower  inlet  velocity  and/or  a  lower  inlet  pressure  so  that  the 
fuel  has  a  sufficiently  long  time  to  react  with  the  catalyst  in 
the  gas-diffusion  layer  at  the  bottom.  With  the  same  reason, 
a  channel  with  a  larger  height  as  well  as  a  gas-diffusion  layer 
with  a  smaller  porosity  will  also  help  enhance  the  fuel  effi¬ 
ciency  (or  shorten  the  channel  length).  Results  show  further 
that  a  catalyst  of  high  activity  is  always  most  desirable  for 
fuel  cell  design  because  it  gives  a  high  power  density  while 
the  accompanied  efficiency  loss  is  virtually  negligible. 

The  above  conclusions,  nevertheless,  are  made  only  in 
qualitative  senses;  namely,  they  give  only  a  trend  for  the 
design  consideration.  To  obtain  useful  data  for  a  quantita¬ 
tive  design,  experiments  must  be  implemented  to  determine 


relevant  parameters  used  in  the  simplified  model.  Particu¬ 
larly,  the  value  of  §  cannot  be  determined  without  going 
through  a  series  of  experiments  in  which  the  porosity  of  the 
gas-diffusion  layer  is  varied. 
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